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Abstract 

In the present contribution we study a PDE system describing the evolution 
psj ■ of a nematic Hquid crystals flow under kinematic transports for molecules of 

^ . different shapes. More in particular, the evolution of the velocity field u is ruled 

by the Navier-Stokes incompressible system with a stress tensor exhibiting a 






CS ■ special coupling between the transport and the induced terms. The dynamics of 



cn 



the director field d is described by a variation of a parabolic Ginzburg-Landau 



t^^ \ equation with a suitable penalization of the physical constraint |d| = 1. Such 

^■"■^ ■ equation accounts for both the kinematic transport by the flow field and the 

internal relaxation due to the elastic energy. The main aim of this contribution 
is to overcome the lack of a maximum principle for the director equation and 
'rj ■ prove (without any restriction on the data and on the physical constants of 

'V^ ■ the problem) the existence of global in time weak solutions under physically 

5^ . meaningful boundary conditions on d and u . 

Key words: liquid crystals, Navier-Stokes system, existence of weak solutions. 
AMS (MOS) subject classification: 35D30, 35K45, 35Q30, 76A15. 

1 Introduction 

In this paper we consider a hydrodynamical system modeling the flow of nematic liquid 
crystals. Assuming that the material occupies a bounded spatial domain fi C M^ 
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with a smooth boundary F, the system couples the 3D incompressible Navier-Stokes 
equations governing the motion of the velocities with a modified AUen-Cahn equation 
for the director field, that is 

divu = 0, in(0,r) x Q, (1.1) 

dtU + dw{u^u) + Vp = dwT + f, in(0,T)xfi, (1.2) 

dtd + u ■ Vd - ad ■ Vu + {1 - a)d ■ V^u = 7(Ad - VdW{d)), in (0, T) x fi, 

(1.3) 

where 

T = § - A (Vd Vd) - aA(Ad - VdW{d)) ® d+ (1 - a)Ad® (Ad - VdW{d)), (1.4) 

S = /i (Vw + V^w) , (1.5) 

T and S being the Cauchy stress and the Newtonian viscous stress tensors, respec- 
tively. Here u denotes the velocity field of the fiow, d is the director field and stands 
for the averaged macroscopic/continuum molecular orientation in M^, p is a scalar 
function representing the hydrodynamic pressure (including the hydrostatic part and 
the induced elastic part from the orientation field) and / is a given external force. The 
positive constants /i , A and 7 stand for the viscosity, the competition between kinetic 
energy and potential energy, and the microscopic elastic relaxation time (Deborah 
number) for the molecular orientation field, respectively. The function W penalizes 
the deviation of the length |d| from the value 1, which is due to liquid crystal molecules 
being of similar size (cf. [H]). A typical example is a double well potential as, e.g., 
W{d) = (|<ip — 1)^. In general W may be written as a sum of a convex part, and 
a smooth, but possibly non-convex one. Finally, the constant a G [0, 1] is a param- 
eter related to the shape of the liquid crystal molecules. For instance, the spherical, 
rod-like and disc-like liquid crystal molecules correspond to the cases a = ^, 1 and 0, 
respectively (cf., e.g., |3], [7] and [2Tj). 

Concerning the notation, Vd represents the gradient with respect to the variable 
d. VdQVd denotes the 3x3 matrix whose {i,j)-th entry is given by VidVjd, for 
^ < "^5 j < 3, and ® stands for the usual Kronecker product, i.e., {u ® u)ij := UiUj , 
for i,j = 1,2,3. Finally, V"^ indicates the transpose of the gradient. 

We notice that this system was very successful in describing the coupling be- 
tween the velocity field u and the director field d, especially in the liquid crystals of 
nematic type. 

The hydrodynamics theory of liquid crystals was due to Ericksen and Leslie (cf. 
[3] and [ID])- However, the general Ericksen-Leslie system was so complicated that 
only some special cases of it have been investigated theoretically or numerically in the 
literature. 

In this context, Lin and Liu (cf. [11] and [12]) formulated a simplified version of 
the original model which has been analyzed also by several other authors, see, e.g., [18], 
|19j . [23] . In the simplified model, some meaningful physical terms, like the stretching 
and rotation effects of the director field induced by the straining of the fiuid, are not 
taken into account. 



In a following paper by Coutand and Shkoller [2], the authors considered a 
model in which the stretching term is present and proved a local well-posedness result. 
Here, due to the presence of the stretching term, the total energy balance does not 
hold. To overcome such an inconvenience. Sun and Liu [21] proposed a variant of the 
Lin and Liu model [H] in which not only the stretching term is included in the system, 
but also a suitable new component to the stress tensor is added. 

In our paper, we refer to a slightly more general model derived by Wu, Xu and 
Liu in 123]. 

The main interests in the topic come essentially from two directions. First, 
the previous results in the literature were only obtained in 2D or in 3D, under the 
assumption that the viscosity coefficient fj, in the stress § (cf. (11.51) ) is sufficiently big 
with respect to proper norms of the initial data and with respect to other coefficients 
like A. In the previous contributions [21] and [23] it was claimed that, due to the 
impossibility of proving the boundedness in L°° for d (because the maximum principle 
cannot be applied to the director field equation), the existence of solution was out of 
reach without assuming to have a big viscosity coefficient in the velocity equation. Our 
main result shows that, even if the existence of classical solutions cannot be proved 
without any restriction on the size of the coefficients and the data (as it is for the 
uncoupled 3D Navier-Stokes system), however, it is possible to obtain the existence 
of weak solutions. This is in agreement with the previous contributions in the field of 
incompressible 3D Navier-Stokes equations. 

The main point here is an appropriate choice of the test functions leading to a 
rigorous weak formulation of the system (cf. the following formulas (I3.8H3.10|) ). Let 
us notice that in the recent paper [1] formal computations are performed in order to 
show the existence of weak solutions for such a problem, but no rigorous definition of 
the weak formulation, as well as no proof of existence of such solutions are given. In 
particular, in our manuscript, choosing properly the space of the test functions in the 
weak momentum equation (cf. (13. 9p ). we obtain well-defined weak solutions. This is 
necessary in order to deal with the stretching term in the Cauchy stress T (cf. (II. 4[) ). 
More comments on this point are given in Remark 13.21 

The second novelty of our analysis consists in the fact that, to our knowledge, 
all the previous contributions in the literature (except for [1] where formal results 
are stated in case of Dirichlet, Neumann and periodic boundary conditions) were ob- 
tained assuming periodic boundary conditions on the director field d. However, from 
the applications point of view, the cases of non-homogeneous Dirichlet or Neumann 
boundary conditions look more appropriate (cf., e.g., [15] where it is pointed out that 
the Neumann boundary conditions for d are also suitable for the implementation of a 
numerical scheme). Here we can rigorously deal with all the three types of conditions 
for d: periodic, Dirichlet and Neumann. 

We note that another new aspect of this contribution relies on the techniques 
employed for the proof of existence of solutions. This method is based on the com- 
bination of a Faedo-Galerkin approximation and a regularization procedure, which is 
necessary in order to treat the high order stretching terms in the weak momentum 
equation. Indeed, a non standard but also physically meaningful regularization of the 
momentum equation is obtained by adding to it an r-Laplacian operator acting on 



the velocities, i.e., we add in the stress tensor a term of the type |V'u|''~^Vii. The 
reader can refer to the series of papers on the J.-L. Lions models (cf., e.g. [13] and 
[m Chap. 2, Sec. 5]) or on the Ladyzhenskaya models (cf., e.g., [8] and references 
therein), where iVwl*" is replaced by |Vit + V'^uY . 

Let us mention that the global in time existence of weak solutions can be con- 
sidered as a starting point for the analysis of the long-time behavior of solutions. This 
problem, up to now, has been considered for this system for instance in [23], where 
the convergence of a global strong solution to a single steady state as time tends to 
infinity has been proved in 2D and in 3D for special sets of data, and in [6] , where the 
existence of a smooth global attractor of finite fractal dimension is obtained. However, 
in both these cases the system was endowed with periodic boundary conditions. 

Finally, our results have been recently used in [TT] and in [1] where the authors 
prove, via Lojasiewicz-Simon techniques, the convergence of the trajectories to the 
stationary states of system fll.lH1.3l) - coupled with suitable boundary conditions - 
and the existence of weak solutions for a non- isothermal system with Neumann (for 
d) and complete slip (for u) boundary conditions, respectively. 

Plan of the paper. In the following Section [2] we briefly introduce the modeling 
approach leading to our system and we discuss the choice of the boundary conditions. 
The weak formulation of (ILlHl.Sp is given in Section [3l where the main theorem 
regarding existence of global in time solutions is stated. The proof is given in the two 
remaining Sections H] and O In particular, in Section H] the a priori estimates, from 
which we deduce a rigorous character of the approximated Faedo-Galerkin scheme 
presented in Section [5], are obtained. 

2 Mathematical model 

In this section we briefly derive system 01.lHl.3p from the macroscopic point of view. 
Hence, suppose that the material occupies a bounded spatial domain fi C M'^ , with a 
sufficiently regular boundary. Let u = u(t, x) denote the velocity in the Eulerian ref- 
erence system. Accordingly, the mass conservation is expressed by means of continuity 
equation (i.e. the standard incompressibility constraint) 

div'u = 0, (2.1) 

which is relevant in the context of nematic liquid crystals. 

In the context of hydrodynamics, the basic variable is the flow map (the particle 
trajectory) x{X, t) , where X is the original labelling (the Lagrangian coordinate) of 
the particle, also referred to as the material coordinate, and x is the current (Eulerian) 
coordinate, which is also called the reference coordinate. For a given velocity field 
u{x, t) , the fiow map is defined by the ODE and initial condition 

Xt = u{x{X,t),t), x{X,0)=X. 

In order to incorporate the properties of the material, we need to introduce the defor- 
mation tensor F such that 

¥{X,t)=dxixiX,t)). 



Combining these two equations with the chain rule formula and defining, in Euler 
coordinates, F (x, t) = ¥{X, t) , we obtain the following transport equation for F 

Fi + (it ■ V)F = VwF. 

Without ambiguity, in what follows we do not distinguish between the two symbols F 
and F . In this case the transport of the director field d can be stated as 

d{x{X,t),t) =F(io(X), 

where do is the initial condition. Taking now the full time derivative of both sides we 
get 

^d{x{X,t),t) = ^do{X) = Vu¥do{X) = Vud={dV)u. 

This allows us to deduce that the total transport associated to the orientation field d 
is 

dfd + u ■ Vd — d ■ Vu 

and so we get the following transport equation for d 

dtd + u ■ Vd - ad ■ Vu + {1 - a)d ■ V'^u = -f{Ad - VdW{d)). (2.2) 

In the general case, W may be a penalty function that can be written as a sum of 
a convex (possibly non smooth) part, and a smooth, but possibly non-convex one. 
Equation (12.21) is associated with conservation of angular momentum. The left-hand 
side stands for the kinematic transport by the fiow field, while the right-hand side 
denotes the internal relaxation due to the elastic energy (cf., e.g., [21]). 

Finally, following the lines of [23], i.e., applying the Hamilton's principle to the 
action functional 



Jo Juo -^ 



^\J'~^VJ'do{X)\' + W{J'do{X)) 



JdXdt, 



where f/o is the region occupied by the fiuid at time t and J = det^ , we deduce the 
following conservation of linear momentum equation 

dtu + div (w (g) n) + Vp = div (fx ( Vu + V'^'u) - A (Vd Vd) (2.3) 

- aX{Ad - VdW{d)) ® d + (1 - a)Xd ® (Ad - VdW{d))) + /. 



This relation combines a usual equation describing the fiow of an isotropic fiuid with 
an extra nonlinear coupling term that is anisotropic. The extra term is the induced 
elastic stress from the elastic energy through the transport, which is represented by 
the equation for d. 

On account of the previous analysis, we get system (I1.1H1.3I) . This system 
must be supplemented with a suitable set of boundary conditions: the homogeneous 
Dirichlet boundary condition for the velocity field 

u = 0, on(0,r)xr, (2.4) 



together with the Neumann homogeneous boundary condition for the director field 

^r^d = 0, on(0,T)xr. (2.5) 

Let us note that with similar techniques we can also treat the case of Dirichlet bound- 
ary condition for d 

d\r = h, on(0,T)xr, (2.6) 

assuming the boundary datum h regular enough. In the following we will detail the 
proper modifications in the analysis in order to treat also this case. 

Finally, we can deal also with the case of periodic boundary conditions, but, 
since in this case the computations are easier, we will not enter into full details. 

3 Main results 

For the sake of simplicity, we restrict ourselves to the case 7 = A = 1 . Accordingly, 
our problem 01.1H1.3p . endowed with initial and boundary conditions, reads as follows 

div'u = 0, in (0,T) x fi, (3.1) 

dtu + div (n (g) w) + Vp = div {fi {Vu + V'^n) ) - div (Vd Vd) (3.2) 

- div (a(Ad - VdW{d)) ® d - (1 - a)d ® (Ad - VdW{d))) + /, in (0, T) x n, 

dtd + u ■ Vd - ad ■ Vu + {1 - a)d ■ V^u = (Ad - VdW{d)), in (0, T) x Q, 

(3.3) 

ti(0, ■) = -uo, d(0, ■) = do, in n, (3.4) 

u = 0, on (0,T) X r, (3.5) 

dr,d = 0, on (0,T) x F. (3.6) 

As already mentioned in the previous section we will also consider the same problem 
with Dirichlet boundary conditions for d, that is 

d|r = /i, on(0,T)xF (3.7) 

instead of condition (13. 6p . 

To begin with, we introduce a weak formulation of (I3.1H3.6P and state our 
main result on the existence of global-in-time weak solutions, without any restriction 
imposed on the initial data or on /i. 

3.1 Weak formulation 

In the weak formulation, the momentum equation (13. 2p together with the incompres- 
sibility constraint (13. ip are replaced by a family of integral identities 

u{t,-)-Vip = 0, for a.a. t G (0,T), (3.8) 

Jn 

{dtU, ip) - u®u:V(p+ fi {Vu + V^u) : Vcp = (3.9) 

Jn Jn 



+ / (Vd Vd) :V(p + a / (Ad - VdW{d)) ®d:Vif 
Jn Jn 

-(1 - a) / d ® (Ad - VdW{d)) : V^ + / / ■ ^, for a.a. t G (0, T), 

for any ip e PFo'^(fi; M^) such that divv? = 0. 

Equation (13. 3p holds in the strong sense, thanks to the regularity obtained for 
d. More specifically, we have 

dtd + u-Vd-ad-S/u + {l- a)d ■ V'^u = Ad - S/dW{d), a.e. in (0, T) x Q, 

(3.10) 

dr^d = 0, a.e. on (0,T) X r, (3.11) 

d(0, ■) = do, a.e. infi. (3.12) 

A wea/i; solution is a pair (w, d) satisfying (I3.8H3.12|) . with u{0, ■) = uq, a.e. in Q. 

A weak solution of the Dirichlet problem is a pair (w, d) satisfying (I3.8H3.10|) . 
(I3.12p . with it(0, ■) = Uq, a.e. in ^2, and 

d\r = h, a.e. on (0,T) X r. (3.13) 

3.2 Main existence theorems 

Before formulating the main result of this paper, let us state the list of hypotheses im- 
posed on the constitutive functions. We assume that fi and a are positive coefficients, 
with a e [0,1], and 

W G C^{R^), W>0, (3.14) 

W = Wi + W2 s.t. Wi is convex and W2 G C\M.^), VW2 G C°'^(M^ M^) (3.15) 
f GL2(0,T;iy-^'2(fi;M3)). (3_ig) 

Our first main result reads as follows. 

Theorem 3.1. Let Q G M.^ be a bounded domain of class C^'^ . Assume that hy- 
potheses Ii3.14[m^) are satisfied. Finally, let the initial data be such that 

UoeW^^\n;R^), div uo = in L\n), (3.17) 

do G W^'^n; R^), W{do) G L\n). (3.18) 

Then problem I\3.8[\3.l2\) possesses a global in time weak solution (u, d) belonging 
to the class 

u G L~(0, T; L\n; R^)) n ^^(0, T; Wq'^Q; R^)), (3.19) 

dtu G L\0, T; W-^'^^\n; R^)), (3.20) 

W{d) G L^{0, T; L\n)), VdW{d) G ^^((0, T) x fi; R^), (3.21) 

d G L°°(0, T; W^'^n; R^)) H L^O, T; W^'^Q; R^)) n H\0, T; L^''^{yt; R^)), (3.22) 



and additionally satisfying, for a.a. t G (0, T) , the energy inequality 

I f {\u\' + \Vd\' + 2W{d)){t) + 2\\{-Ad + VM{d)){t)\\l,^^,^^,^ (3.23) 

+ /ill Vw(t) 11^2(^.^3X3) < C'||/(t)||^-l,2(f^.]R3), 

where C denotes a positive constant depending on Q . 

Remark 3.2. The regularity of the test function ip in (13.91) can be justified by noting 
that, thanks to (13191) . (l3:2T|l and (131221) . we have 

u G L°°(0, T; L^iVl- M^)) n ^^(0, T; W^'^iVl- M^)), 

Vd G L°^(0, T; ^^(f]; M^^^)) n ^^(0, T; Vr^'2(fi; R^^^)), 
Ad - WdW{d) G L2((0, T) X fi; M^)^ de L°°(0, T; ^^^'^(fi; M^)), 
and hence 

u®u, VdQVd, {Ad-VdW{d))®deL\0,T;L^/'^{n;R^'''^)), (3.24) 

whence their (distributional) divergence are elements of the space 

The same problem (cf. also [21] and [23]) occurs in the 2D case, where the same weak 
formulation is needed in order to get existence of weak solutions. 

Finally, we can state the second result in case of Dirichlet boundary conditions 
for d. 

Theorem 3.3. Let Q cM^ be a bounded domain of class C^'^ . Assume that hypothe- 
ses I13.14[WT5\) and I13.17[WW\) are satisfied. Suppose moreover that the assumptions 

heH\0,T;H-^/\T;M:^))nL°^{0,T;H^/\T;R^)), h{0) = do\r (3.25) 

hold true. Then problem I13.8[]3.1(J\] . Il3.12[m^) possesses a global in time weak solu- 
tion (u, d) belonging to the class stated in / l^.l^ff^T^) . 

The rest of the paper is devoted to the proofs of Theorems 13.11 and 13.31 
In Section m we will prove the a-priori bounds on the solutions, detailing in Re- 
mark 14.11 the differences between the Neumann and the Dirichlet cases. In Section \5\ 
we introduce the regularization schemes and we perform the estimates on the approx- 
imated solutions in order to pass to the limit in the approximated solutions. Since 
in the case of Dirichlet boundary conditions (I3.13P (instead of (13.111) ) the technique 
is analogue, in Section \5\ we perform the estimates and the approximation-passage to 
the limit procedure on the Neumann system (I3.8H3.12T) . 
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4 A priori bounds 

We establish here a number of formal a priori estimates. These will assume a rigorous 
character in the framework of the approximation scheme presented in Section [5] below. 
We inform the reader that a similar technique has been used also in the subsequent 
work [1] (with respect to the present one) for a non-isothermal model with different 
kind of boundary conditions. 

Take (p = u in (jM]) and test fl330|l by -Ad+VdW{d) on Q. Summing up the 
two resulting equalities, using the divergence theorem together with (13. 8|) . we obtain 



15 /■/, |2 



Vd\^ + 2W{d))+fi f \Vu\^+ f \-Ad+VdW{d)\^=H-i {f,u)^i,2. 
Q Jn Jn ° 

(4.1) 

Moreover, applying Schwarz and Poincare inequalities on the right hand side, we can 

deduce the energy estimate (13.231) . Integrating over (0,T) the inequality (14. ip . and 

using assumption (I3.16p . we get the a priori bounds 

u e L°°(0, T; L^(n; M^)) n /.^(O, T; W^'^i^; R^)) n L^°/^((0, T) x fi; M^), (4.2) 

deL°°{0,T;W^'^{n;R^)), (4.3) 

- Ad + VdW{d) G L2(0, T; ^^(fi; M^)) _ (4.4) 

From (14. 4p . on account of (12. 5p and (I3.14p . we get 

/ [\Ad\^+f f V{VdW{d))Vd= f [mAd, 
Jo Jn Jo Jn Jo Jn 

ra being a function in L2(0, T; L'^{n; M?)) . 

Using once more assumption (13.151) . from the previous equation we get 



\ r I |Adp < r I \V{VdW,{d)md\ + \r I imp. 
■^ Jo Jn Jo Jn ^ Jo Jn 



'0 Jn Jo Jn ^ Jo Jn 

Recalling that |Vd| G L°^{0,T; L'^{Q;R^)) (cf. (gSD), then it holds (see again (^^ 
and (Km ) 

deL'^{0,T;W'^'\n;R^)), VdW{d) e L\{0,T) x n;R:^). (4.5) 

From this result, since it-VcZ and dVu belong to L^(0, T; L^'^(r2; R?)) , by comparison 
with (I3.10p we deduce 

dtdeL'iQ.T-L^'^n-W^)). (4.6) 



Now, choosing g(l — a) = 2 in the following interpolation inequality 



^"llL»(f7;IR3x3) ^ t^lll Va||^2(^lg3x3)|| VCi||^6(QK3x3) , l^-'J 



holding true for 

s, gG [l,+oo), aG (0,1), 1/s = (1 - a)/6 + a/2 , (4.8) 
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and using ( l4.3H4.5p . we get 

which, taking s = 10/3 , gives in particular 

Vd G L^°/=^(0, T; Li°/3(fi; M=^^=^)). (4.10) 

This estimate turns out to be crucial for the proof of existence of solutions. Hence, as 
a consequence of the previous estimates, we get 

(- (Vd Vd) + a{Ad - VdW{d)) ® d - (1 - a)d ® (Ad - VdW{d))) (4.11) 
eL5/^((0,T) xfi;M3x3) 

and 

(- (Vd Vd) + a(Ad - VdW{d)) d - (1 - a)d (Ad - VdW{d))) (4.12) 
gL2(0,T;L3/2((];M3x3)). 

Observe that from the a priori estimates (14.21) . (14.31) and (14.61) we derive the solu- 
tion regularity classes (I3.19P and (I3.22p . from which it follows (I3.24p and then (13.200 . 
Moreover, we can prove that the set of the (weak) solutions to problem (I3.8H3.101) is 
weakly stable (compact) with respect to these bounds, namely, taking any sequence 
of (weak) solutions satisfying the above uniform bounds then it admits a convergent 
subsequence. We omit the proof of the weak sequential stability, leaving the details 
to the reader, and we devote the following section to the proof of Theorem 13.11 More 
precisely, we will construct a suitable family of approximate problems whose solutions 
weakly converge (up to subsequences) to limit functions which solve the problem in 
the sense specified in Subsection 13. 1[ 

Remark 4.1. Let us detail here the different energy estimate we obtain in the Dirichlet 
case. If we we take non- homogeneous Dirichlet boundary conditions for d, that is 
fl37[3|) instead of fl3ll]) . we get 



ld_ 
2di 



I (|up + |Vd|2 + 2W^(d)) + /i /" |Vn|2 + f \-Ad + VdW{d)\^ (4.13) 
Jn Jn Jn 



=/f-i/2(r) {ht,dnd)fji/2(^r)+H--^ {fi'^)w^-^- 



Assuming h satisfying (I3.25p . we can handle the first term on the right-hand side of 
(I4.13P using standard trace theorems and regularity results for elliptic equations (cf., 
e-g-! [ISl Lemma 3.2, p. 263]) in this way 

H-^/'2(T){ht,dnd) ^i/2(Y'^ < C*||^t||i:/-i/2(r)||d||//2(f7) (4-14) 

<r(\\h\\'^ -^ll/ll|2 W"'"IIAWI|2 

^ '-' \\\'^t\\H-^/-2{T) + ll"'llii'3/2(r)y + Tll^"llL2(n) • 
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In order to treat the last term in fl4.14p we estimate from below the last integral in 
IKT^ as follows 

II - Ad + VdW^(d)||i.(^) = ||Ad||i.(^) + ||Vdl^(d)||i.(^) -2(Ad, V^W-(d)) (4.15) 

>||Ad||i,(^. + ||VrfW^(d)||i.,^. 



+ 2 [ V{VdW{d))Wd-2 f dnd{VdWid))ir 
Jn Jv 

>||AcZ||i.(^)-C||Vd||i.(^)-i||Ad||i.(^) 



4' 



C||^|lWrn-C||V.iy(/i)||i2 



Here we have employed assumptions fl3.14p . (13.151) and again standard elliptic esti- 
mates and trace theorems. Integrating over (0,T) the inequality (I4.13p . and using 
assumptions (I3.14p . (I3.16p . (I3.25P together with (I4.14H¥]T^ and a standard Gronwall 
lemma, we get the a priori bounds 

u e L°°(0, T; L^{VL] M^)) n L2(0, T; W^^'^{VL] M^)) n L^°/^((0, T) X Vt] M^), 

-Ad + VdW{d) e L2(0, T; L\n; M^)), 

d G L\0, T- H\n; M^)), VdW{d) G ^^((0, T) x (]; M.^) . 

Note that we have used, in particular, the assumption h G L°°(0,T; if^/^(r; M'^)) in 
order to get h G L°°(0,T; C°(r;M3)) , which, thanks to the fact that W G ^^(M), 
implies \\V dW{h) Wl^^^^ G 1/^(0, T). The remaining bounds dHUHHSD still holds true 
in the Dirichlet case. 



5 Approximations 

Here we introduce a double approximation scheme: a standard Faedo-Galerkin method 
is coupled with an approximation of the convective term and a regularization of the 
momentum equation by adding an r-Laplacian operator acting on the velocities. Re- 
garding the choice of these two regularization, let us note that for the first one (related 
to the convective term) we follow the classical approach by Leray [n], while the mo- 
mentum equation with the r-Laplacian acting on the velocity field u (cf. the following 
(15.11) ) is a meaningful approximation of the standard Navier-Stokes system. Indeed, 
one can refer to the series of papers on the J.-L. Lions models [13] and [T^ Chap. 2, 
Sec. 5] or on the Ladyzhenskaya model (cf., e.g., [8] and references therein), where 
|Vw|^ is replaced by \Vu + V'^u\^ . 

The standard technique of the Faedo-Galerkin approximation scheme will be ap- 
plied to construct the solutions to the Navier-Stokes system (I3.8ti3.9p (see Temam [22]). 
To this aim, let us introduce first the Hilbert space 

W^^^^ = {ve Wo'^(fi;R2) I divt; = 0, a.e. in fi} 
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and take an orthonormal basis {Vn}'^=i ■ We fix M, N & N such that M < N and 
consider the finite-dimensional space X^ = spa.n{Vn}n=i- 

Then the approximate velocity field Un,m G C*^([0,T]; Xjv) solves the Faedo- 
Galerkin system 

-77 / Un,M V = [un^m]m ® Unm • '^^ ~ TJ / |V'Uiv,Mr"^V'UAr,A/ ■ Vv (5.1) 



fllVuN, 



M + V^un,ai ] ■.Vv+ I VdN,M VdNM : Vv 



n 



+a / (Adyv,M - VdW{dN,M)) ® djv,A/ : Vt; 



n 



-(1 - a) / d;v,A/ ® (AcZ^,A/ - VdW{dN,M)) ■.Vv+ [ fv, for all t E [0, T], 



UN,Mi^i ■)-v = / Wo ■ f , (5.2) 

Jn 

for any v G Xjy and r G (3,10/3). Here, the symbol [v]m denotes the orthogonal 
projection onto the space Xm = spaii{vn}n=i- We observe that we need to intro- 
duce the additional term jj\VuN,M\^~'^'Vui\f^M (cf. (13.91) ) in order to obtain enough 
regularity for the velocity field in the director equation. 

In (15. ip appears also the function d^^M which is determined in terms of un^m 
as the unique solution of the system 

dtdN,M + un,m ■ ydN,M - adNM ■ yuN,M + (1 - a)dN,M ■ V^un,m + VdW{dN,M) (5.3) 

= AdN,M, in(0,r)xf], 

5r^d^,A/ = 0, on(0,T)xr, (5.4) 

<i7v,Af (0, ■) = cio,A//, in fi, (5.5) 

do,M being a suitable smooth approximation of the initial datum do (cf. (I3.10HHTT^ ). 

We will follow the original approach to the Navier-Stokes system by Leray [S] in 
order to regularize the convective terms in (15.11) . (15.31) . Hence, for any fixed M, N , we 
can solve problem (I5.1H5.5I) by means of a fixed point argument, exactly as detailed 
in [5l Chapter 3]. 

Indeed, observe that all the a priori bounds derived formally in Section H] still 
hold for our approximate problem. Hence, if we fix tt G C{[0,T];Xn), then we can 
find d = d[u] solving (I5.3H5.5I) . Inserting d[u\ in system (I5.1H5.2I) we can define a 
mapping u i— )■ T[u\ , T[u\ being the solution of the system. On account of the a priori 
bounds obtained in Section HI we can easily show that T admits a fixed point by means 
of the classical Schauder's argument on (0,To), with < Tq < T. Finally, applying 
again the a priori estimates, we are allowed to conclude that the approximate solutions 
can be extended to the whole time interval [0,T] (see [H Chapter 6] for details). 

Now our strategy consists in passing to the limit first for TV — )■ oo and then for 
M -^ oo. This will be explained in the following subsections. 
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5.1 Passage to the limit as A/^ — )> oo 

On account of the regularizing term introduced in fIS.ip . from the corresponding energy 
estimate we now obtain 

M'i Viijv,MirL.((o,T)xn;R3x3) < C, for T G (3, 10/3), (5.6) 

from which we infer that, for any fixed M, the set of functions \Wum^mY~'^'^Un,m is 
uniformly bounded in L^((0,T) x fi;]R3^2). Note that, since r G (3, 10/3), it holds 
r/(r-l)G (10/7,3/2). 

Hence, we deduce the following convergence results 

Unm -^ Um weakly-(*) in L°°(0, T; L^iVt; M^)) n ^^(0, T; W^''^{yt- M^)), (5.7) 

^Unm -^ Vwa/ weakly in ^'^(O, T; U{VL] M^)), (5.8) 

dtUN,M -^ dtUM weakly in L^(0,T; W'^'^'/'-^n; R^)), (5.9) 

dN,M -^ du weakly-(*) in L~(0, T; W^''^{yi- M^)) n ^^(0, T; ^^^'^(fi; W^)). (5.10) 

Observe that in (15. ip the projection on Xm is kept in the convective term when passing 
to the limit as A^ — )■ oo . 

Moreover, by virtue of (I5.10p and a simple interpolation argument, we have also 

Vdjv.A/ -^ VdM strongly in L^((0, T) x Vt; M^^^), for r/ G [1, 10/3). (5.11) 

Going back to (15. 7p and (15.81) . by standard interpolation results, some embedding 
properties of Sobolev spaces, and the Aubin-Lions lemma, then we get 

Un,m ^ Um strongly in L''((0,T) X 1];]R^), for some s > 5. (5.12) 

A combination of (15.121) with (15. lip gives 

un,m ■ VdAT^A/ — ^ Um ■ ydM strongly in L''((0,T) x f2), for some s > 2, (5.13) 

whereas, combining (15. 8 p and (I5.10p . one obtains 

dN,M ■ ^un^m — ^ dM • ^um weakly in L^((0, T) x (]), for some p > 2. 

Moreover, it holds 

dtdN,M -^ dtdM weakly in L2(0, T; L^{Vl- M^)), (5.14) 

and we have also 



|Vn^,Afr^^Vn7v,M -^ \VumY~'^Vum weakly in L''/"-i((0,T) x n-.W^""^). 

Considering now the pair of limit functions {uM,dM), it can be proved that it solves 
the problem 

/ WM(i, ■) ■ Vy? = 0, for a.a. tG (0,T), (5.15) 
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/ {dtUM,'f) - ([um]m ® Um : Vv?) + / / yu (Vum + V^Um) : Vy? (5.16) 

JO Jo Jn ^ ' Jo Jn 

= (Vdjv/ VdM + a (AcZm - VdW^(cZM)) ® ^a/) : Vy^ 

- f / (1 - a)cZM ® (AcZm - VdW{dM)) : V(^ 
JO Jn 

-i- / /|VuM|'-'Vnj;,-Vv9+ / //■(/^, foralHG(0,T), 
^^^ Jo Jf7 Jo Jn 

for any ^ G C°°{n; R^) such that div(^ = 0. 

Passing to the hmit as A^ — ?■ oo also in the equation for d, we get, a.e. in 

(o,T) X n, 

dtdM + Um ■ VdM - adM ■ Vum + (1 - a)dM ■ V^um = ^dM - VdW{dM), (5.17) 

as well as 

dr,dM = 0, a.e. in (0,T) X r, (5.18) 

<iA/(0, ■) = cZo,A//, a-e. in Q. (5.19) 

Going back to (15. 11) . if we replace v by ■uat^a/ and then integrate with respect to time 
over (0,t), we obtain 

\\uN,Mi't)\\h(n)+ / / yt^|VWAr,A/ + V^ttAT^Mp + TJ / / \'^Un,m\'' (5.20) 

Jo Jn ^^^ Jo Jn 

= ll^o|li2(n) + 2 / / {VdN,M VdN,M) : Vwat^m 
Jo Jn 

+2a / / {AdN,M - Vdiy(div,M)) ® ^at^m : Vhat^m 
Jo Jn 

-2(1 - a) / / dN,M (AcZat^m - VdVr(djv,M)) : Vwat^m, + / / / • -uat^m, 
Jo Jn Jo Jn 



for all t G (0, T) . 

Next, on account of fl5.7H5.9l) . we can consider Um as a test function in fl5.16p 
and get 

W^Mimhin) + j j /^iVi^A/ + VV/P + Jfjj l^^^r (5.21) 

= ll^olli2(n) + 2 / / (AcZm - VdVr(dA/)) ^M : Vum 
Jo Jn 



+2a / / {AdM - VdW^(dA/)) dM : Vi^m 
Jo Jn 

-2(1 - a) [ j dM®{AdM-VdW{dM)):VuM + 2 f [fuM, for alH G (0,T). 
Jo Jn Jo Jn 
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Actually, the L'' -regularity (15.81) of Vum is essential at this level since the terms 
{AcIm — VdVr(dM)) 8' dM and (Ad^/ — "^dWidiM)) ® d^ do not necessarily belong 
to L^ . More precisely, the best we can conclude it is that they lie in L^'^ (cf. fIS.lOp 
and fIS.lip ). Now, multiplying fl5.3p by Ad^^M — VdW^d^^M) and integrating on 
(0, t) X fl, we obtain 

II Vd;v,A/(t)||i2(^) + 2 / Wid^^Mm + 2 [ f \Ad^,M - Vrfiy(d^,M)r (5.22) 
Jn Jo Jn 

= II Vdo,M|li2(n) + 2 / W{do^M) + 2 / {un,m ■ Vd^,A/, AdN,M - ^dW{dN,M)) 
Jn Jo 

+2 / [-adN,M ■ Vun,m + (1 - a)dN,M ■ V^un,ai, Ad^M - VdW{dN,M)) , 

for all t G (0,T). If we multiply (15.17P by AcZm — VdW{dM) and we integrate on 
(0,t) X fi, then we deduce (for any t G (0,T)) 

l|VdM(t)||i2(f,) + 2 /" W^(dM)(t) + 2 / / |AdM - VdW{dM)\^ (5.23) 

in JO Jn 

= ||Vdo,M||i2(f,)+2 f W{do,M) 

Jn 

+2 / (um ■ ydM - adM ■ Vum + (1 - a)dM ■ V^Um, Ad^/ - WdW{dM)) ■ 
Jo 

Note that this is possible since fl5.17p makes sense as a relation in L^ (cf. (15. 1 SI - 
IS. 14p ). thanks to the higher regularity (15. 8p and (15.121) of um and Vu^v/ given by the 
regularizing term ■jj\WumY~'^'^Um in (15. ip . Summing (I5.20p with (I5.22p and then 
(I5.2ip with (15. 23 p . passing to the limit as A^ ^^ oo we obtain 

/ / |VwAr,A/r -^ IVwa/T'^^Vum : VtiAf, 

Jo Jn Jo Jn 

[ [ \Adr,,M - VdW{d^,M)\' ^ [ [ l^dM - VdW{dM)\'. 

Jo Jn Jo Jn 

So that, by means of standard Minty's trick and monotonicity argument, we infer 

Vun,m -> Vwa/ strongly in L'^((0,T) x fi;M3''3), 

AdN,M -^ AdM strongly in L^((0,T) x Q; R^). 
This concludes the passage to the limit as A^ — )> oo . 
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5.2 Passage to the limit as M — )> oo 

The final step in the proof of the main result consists in passing to the limit as M — )■ oo 

in dSlsHSISD. 

First, we observe that we can still deduce the convergence results in (15. 7p and 
(I5.12p when taking M — )■ oo . Moreover, the following convergence results hold true 

dtUM -^ dtu weakly in L^(0,T; iy-^''^/"-^(l]; M^)) , (5.24) 

dM -^ d weakly-(*) in L~(0, T; W^^'^iVt- M?)) n L2(0, T; W'^^'^iVt- M?)) , (5.25) 

dtdu -^ dtd weakly in L2(0, T; L^'^iVt- M^)) , (5.26) 

and in particular 

M-^/^-'-^^Vum -^ strongly in ^'-^((O, T) x Vt- M^^^^). (5.27) 

We are now in a position to pass to the limit as M — > oo in (l5.15f[5?T9|) and finally 
recover ([SiHHSll]). 

Acknowledgements. We would like to thank Sergio Frigeri, Gianni Gilardi and the 
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